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Abstract 

A number of approaches to the arrival time problem employ a complex potential of a simple 
step function type and the arrival time distribution may then be calculated using the stationary 
scattering wave functions. Here, it is shown that in the Zeno limit (in which the potential becomes 
very large), the arrival time distribution may be obtained in a clear and simple way using a path 
integral representation of the propagator together with the path decomposition expansion (in which 
the propagator is factored across a surface of constant time). This method also shows that the 
same result is obtained for a wide class of complex potentials. 
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I. INTRODUCTION 



Some of the interesting outstanding problems in quantum theory concern situations in 
which time appears in a non-trivial way. One of the simplest such problems is the arrival 
time problem. In the one-dimensional version of this problem one considers an initial wave 
function concentrated in the region a; > and consisting entirely of negative momenta. The 
question is then to find the probability that the particle crosses a; = between time r and 
T -|- dr. The classical analysis of this problem is trivial, but the quantum analysis is not so, 
due in part to the fact that the usual machinery of quantum measurement refers to fixed 
moments of time and not to measurements distributed over an interval of time. 

There are many different approaches to this problem [1]. One particular approach is to 
include a complex potential 

V{x) = -iVoe{-x) (1) 
in the Schrodinger equations, and to then compute the final state 

mr)) = exp (^-^HoT - ^0{-x)t^ (2) 

where Hq is the free Hamiltonian. The intuitive idea behind this is that the part of the wave 
packet that reaches the origin during the time interval [0, r] is absorbed, so that 

7V(r) = (^(T)|^(r)) (3) 

is the probability of not crossing during the time interval. The probability of crossing 
between r and t + dr is then 

n(r) = -- (4) 

Such potentials were originally considered by AUcock in his seminal work [2] and have 
subsequently appeared in detector models of arrival times [3, 4]. A recent interesting result 
of Echanobe et al. is that under certain conditions a complex potential of the form Eq.(l) 
is essentially the same as pulsed measurements, in which the wave function is measured at 
discrete time intervals [5]. 

The difficulty behind this approach, however, is that the wave function is not entirely 
absorbed by this complex potential and N{t) includes parts of the initial wave function that 
have refiected off the potential. The refiection is small for small Vq but then the resolution 
of the measurement, which is proportional to h/Vo, is poor. On the other hand, there is a 



lot of reflection for large Vq and indeed the wave function is entirely reflected in the limit 
Vo — > oo. This is the quantum Zeno effect [6] and plagues many different approaches to the 
arrival time problem. 

Echanobc et al. have made an interesting proposal which embraces the Zeno effect in 
the Vq — * oo limit yet at the same time extracts the physics hidden within it by suitable 
normalization [5]. They consider the limit Vq ^ oo of the expression 

which is normalized since dr n(T) = 1 — N{oo). The point is that N{oo) represents the 
total amount of reflected wave function, so 1 — N{oo) represents the total probabihty of 
crossing during the time interval [0, oo) and this goes to zero as Vq oo. The expression 
n(r) also goes to zero as Vq — > cxd but the ratio Eq.(5) is finite and independent of Vq, and 
defines a reasonable normalized arrival time distribution function. Known results on the 
stationary scattering wave functions [7] yield the result 

n(^) = ^,/^yi/2 (Mr)mmMr)) (6) 
from which the normalized result is 

^N{r) ^ :^{Mr)mx)p\Mr)) (7) 

where |'«/'/(t)) is the freely evolved wave function and (p) is the average momentum in the 
initial wave packet [5, 7]. 

An interesting question in these expressions is the origin of the form of the p5{x)p term, 
which is not obvious from the calculation of it in Ref . [7] . Compare this to the simplest guess 
for the arrival time distribution function, the current density 

"^^^^ ^ 2^^^^^^^! ^^^^^^ + ^^^^^^ 1^^^^^^ 
which is sensible classically but is not always positive in the quantum case [8]. (Here d{x) = 
|0)(0|, where |0) denotes the eigenstate \x) of the position operator at a; = 0.) A simple 
operator re-ordering of J{t) gives the "ideal" arrival time distribution of Kijowski 

^K{r) = ^{Mr)\p'^'mp'^'\Mr)) (9) 

fiL 

which is clearly positive, but harder to relate to specific measurement schemes [9]. It would 
be of value to find a derivation of Eq.(7) which gives some insight into the form of this 



expression. A further question concerns the role of the complex potential. Although Eq.(7) 
is independent of Vq, it is not clear to what extent the result depends on the specific choice 
of complex potential, Eq.(l). 

In this paper, we show that both of these questions are answered very simply using 
path integral methods. We show that the general form (6) is derived very easily (up to a 
constant, fixed by normalization). Secondly, the result (7) is seen to be true for a wide class 
of potentials of the form 

V{x) = -iVo0{-x)f{x) (10) 

where /(x) is any positive function. 

With the general complex potential Eq.(lO), the arrival time distribution (4) is given by 

n(r) = 2{^^\V{x)\^r) 



= 2-^ y dx f{x)\i^{x,T)\' (11) 
where ip{x,T) is defined in Eq.(2). The key is therefore to evaluate the propagator 

g{xi,T\xo,0) = (xilexp (^-^Hqt - y^6'(-£)/(x)r^ \xo) (12) 
for xi < and Xq > 0. This may be calculated using a sum over paths, 

g{xi,T\xo,0) ^ J'Dxexp(^^S^ (13) 

where 

dt(^^mx^ + iVo0i-x)f{x)^ (14) 

and the sum is over all paths x{t) from x{0) = Xq to x{t) = Xi. 

To deal with the step function form of the potential we need to split off the sections of the 
paths lying entirely in x > or x < 0. The way to do this is to use the path decomposition 
expansion or PDX [10-12]. Each path from > to xi < will typically cross x — 
many times, but all paths have a first crossing, at time ti, say. As a consequence of this, it 
is possible to derive the formula, 

ih r dgr 
g{xi,T\xo,0) = — J dti g{xi,T\0,ti)-^{x,ti\xo,0)\^^^ (15) 

Here, gr{x,t\xQ,0) is the restricted propagator given by a sum over paths of the form (13) 
but with all paths restricted to x{t) > 0. It vanishes when either end point is the origin but 



its derivative at a: = is non-zero (and in fact the derivative of Qr corresponds to a sum 
over all paths in x > which end on x = [12]). It is also useful to record a PDX formula 
involving the last crossing time ^2, 

g{xi,T\xo,0) ^ -— J dt2 -^{xi,T\x,t2)\^^Q g{0,t2\xo,ti) (16) 

These two formulae may be combined to give a first and last crossing version of the PDX, 

^(xi,t|xo,0) = ^J dt2 J dti -^{xi,T\x,t2)\^^^ g{0,t2\0,ti) -^{x,ti\xo,0)\^^^ 

(17) 

The restricted propagators in the two regions arc easier to work with than Eq.(12) since 
the potential is zero throughout x > and —iVQf{x) throughout a; < so these may be 
calculated by standard methods, without the complication of the 9{—x) term. The problem 
of calculating the propagator Eq.(12) therefore essentially reduces to the easier problem of 
calculating it between two points lying on x — 0. This can sometimes be evaluated by 
a mode sum calculation [13], even though the full propagator Eq.(12) is not necessarily 
calculable in this way. 

Returning to the first crossing PDX, Eq.(15), gr is the restricted propagator for the free 
particle, which is given by the method of images expression 

gr{xi,T\xo, 0) = 9{xi)9{xo) {gf{xi, t\xo, 0) - gf{-xi, t\xo, 0)) (18) 
where gf denotes the free particle propagator. It follows that 

^{x,ti\xo,0)l^, = 2||(0,ti|xo,0)^(xo) (19) 
Inserting this in Eq.(15), and rewriting it as an operator expression, we obtain the result 
{xi\ exp (^-^Hqt -^e{-x)f{x)T^ \xo) 
= - ^J\t{x,\ex^p(^-^Ho{T-t)-^9{-x)f{x){T-t)^ 

X S{x)p exp (^-^Ho?j \xo) (20) 

We note the appearance of the combination 5{x)p which clearly corresponds to a "crossing 

operator" (and in derivations of the PDX arises directly from a derivative of 9{x) [11]). Now 
note that the operator d{x) has the simple property that for any operator A 

S{x)AS{x) = S{x){0\A\0) (21) 



This property together with Eq.(20) inserted in Eq.(ll) yields 



n(r) 




(22) 



where H = Ho — iVo9{—x)f(x) is the total (non-hermitian) Hamiltonian. Now the key point 
is that as Vq — >• oo, the integrals over t and f are strongly concentrated around r, so we 
obtain 



Furthermore, by the change of variables s — t — t, s' — t — t' , it is easily seen that 
C is independent of r for large Vq, so is just a constant. It is also easily seen that the 
approximation (23) will hold for Vq ^ k^/2'm, where is the largest momentum in the 
initial state. 

Eq.(23) is the main result and is of precisely the form Eq.(6), up to an overall constant. 
Since Eq.(7) is obtained by normalization of Eq.(6), wc therefore obtain Eq.(7). It holds for a 
general class of potentials, of the form Eq.(lO), not just the simple case f{x) — 1 considered 
in Refs.[5, 7]. The reason for this is that for large Vo the paths in the path integral are kept 
out of the region x < 0, only entering it at just before the final time, so the result has very 
hmited dependence on the detailed form of the potential in x < 0. The result Eq.(23), and 
in particular the p6{x)p form, follow as an almost immediate consequence of the PDX for 
large Vq, together with the simple property Eq.(21). 

To fully verify the validity of this path integral approach, we now calculate the constant 
C in the case f{x) — 1 and look for detailed agreement with the unnormalized result Eq.(6). 
Eq.(24) may be written. 




(23) 



where 




(24) 




(25) 
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where, after a change of variables s — t — t, 

(f){x) ^ j^ds {x\ exp [^^^Hqs - ^e{-x)s^ |0) (26) 

The integrand may be represented as a sum over paths from x = at time s = to the 
final point x < at time r. We may use the last crossing PDX, Eq.(16), to split this into a 
sum over paths from a; = at s = 0, to the last crossing ata; = 0, s = and from there 
propagating entirely in x < to the final x at time s — r. This reads 

= ds dii (x|exp ^-■^i/o(s-ii) - y(s-ii)^p|0)(0|exp ^i/M^ |0) 

(27) 

In the second bracket expression, H is the total Hamiltonian so still includes the 6{—x) 
potential. This is the propagator along the edge of an imaginary step potential which, 
fortunately was calculated in Ref. [13] using a mode sum method (for the case of a real step 
potential, but this is readily continued to imaginary values), and we have 

(0|exp(->)|0)^(5^)"^il^;5^ (28) 

Noting that 

ds du= du ds (29) 

JO Jo Jo Ju 

and changing variables to v — s — u, we see that for large Vq the dominant contribution 
comes from close to = and v — 0. We may therefore let r — > oo, and the integrals then 
factor into a product 



[X = 

m 



f*' (^i-p(4i^o.-^.)pio> f * {£^y" ii^^ (30, 

The first integral may be evaluated using the familiar formula [14] 

with E — iVo, and then applying —ihd/dx. The second is evaluated using the formula. 



^0 



I dx^—^^2V^ (32) 



We thus obtain 

1/2 



—j exp(^ ^VmVo\x\j 



(/>(x) = — exp — y/mVo\x\] (33) 



7 



Inserting in Eq.(25), this gives the final result 

2 

The path integral method used here therefore gives precise agreement with the earlier result 
Eq.(6) obtained by scattering methods. 
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